* WAVE 1 Panama
* 1999-2004
clear all
cd "/Users/chuangchen/Library/CloudStorage/OneDrive-UniversityofPittsburgh/Pela project/PELA DATA/Panama/data_analysis/"
import spss using "BASEDATOS_PANAMA_41.sav"
keep nestu pais partido legis p58 p59 p60 p60a p62 p63 p67 p65 p28 p64 p34 p2901 p2902 p2903 p2904 p2905 p2906 p2907 p2908 p2909 p2910

elabel variable (*) ("")

label variable nestu "Study number"
label variable pais "Country name"
label variable partido "Party name"
label variable legis "Legislature"


label drop labels0
label values pais
tostring pais, replace
replace pais="Panama"

replace legis=2004

gen party = "."
replace party = "PN_PRD" if partido ==1
replace party = "PN_PA" if partido ==2
replace party = "PN_Other" if partido ==3

drop partido 

rename party partido

label variable p58 "Ideology"
rename p58 ID1
replace ID1 = . if ID1 > 10
label variable p59 "Ideology of your party"
rename p59 ID2
replace ID2 = . if ID2 > 10

label variable p62 "Female (dummy 1 if woman)"
rename p62 female
recode female (1=0)(2=1)
label define female_label 1 "female" 0 "male"
label value female female_label

rename p60 religious
recode religious (2=0)
label variable religious "Religious (dummy 1 if believer)"
label define religious_label 1 "believer" 0 "nonbeliever"
label value religious religious_label

gen rel_evangelical=1 if p60a==3
replace rel_evangelical=0 if rel_evangelical!=1
label variable rel_evangelical "Evangelical (dummy 1 if evangelical)"
label define rel_evan_lab 1 "Evangelical" 0 "Other"
label value rel_evangelical rel_evan_lab

gen rel_catholic=1 if p60a==1
replace rel_catholic=0 if rel_catholic!=1
label variable rel_catholic "Catholic (dummy 1 if catholic)"
label define rel_cath_lab 1 "Catholic" 0 "Other"
label value rel_catholic rel_cath_lab

gen rel_other=1 if p60a!=1 & p60a!=3
replace rel_other=0 if rel_other!=1
label variable rel_other "Religious Other (dummy 1 if other)"
label define rel_other_lab 1 "Other" 0 "Evangelical/Catholic"
label value rel_other rel_other_lab

drop p60a

* no question about church attendance
* rename p60b religiosity

label variable p63 "Age (years)"
rename p63 age
replace age = . if age == 99

rename p67 education
label variable education "Highest level of education"
replace education = . if education > 6

rename p65 val_abortion
label variable val_abortion "Opinion about abortion"
replace val_abortion = . if val_abortion > 10

rename p64 val_divorce
label variable val_divorce "Opinion about divorce"
replace val_divorce = . if val_divorce > 10

rename p28 eco_regulated
label variable eco_regulated "Regulate Economy"
label drop labels32
replace eco_regulated = . if eco_regulated > 5
replace eco_regulated = eco_regulated*2
* attention to the rescaling

rename p34 eco_tax
label variable eco_tax "Whether legislator prefers direct or indirect taxes"
label define tax_lab 1 "direct" 2 "indirect"
label value eco_tax tax_lab
replace eco_tax = . if eco_tax !=1 & eco_tax !=2                                                             
rename p2901 eco_prices
label variable eco_prices "Price control"

rename p2902 eco_education_p
label variable eco_education_p "Free primary education"
* combines primary and secondary education

rename p2907 eco_education_u 
label variable eco_education_u "Free university education"
recode eco_education_u (1=1) (2=3) (3=5) (4=7)
replace eco_education_u = . if eco_education_u > 7

rename p2903 eco_housing
label variable eco_housing "Subsidize housing"

rename p2904 eco_employment
label variable eco_employment "Employment creation"
recode eco_employment (1=1) (2=3) (3=5) (4=7)
replace eco_employment = . if eco_employment > 7

rename p2906 eco_health
label variable eco_health "Whether the state should provide health services"
recode eco_health (1=1) (2=3) (3=5) (4=7)
replace eco_health = . if eco_health > 7

rename p2908 eco_unemployment
label variable eco_unemployment "Unemployment insurance"
replace eco_unemployment = . if eco_unemployment > 4

rename p2905 eco_pensions
label variable eco_pensions "Pension"
recode eco_pensions (1=1) (2=3) (3=5) (4=7)
replace eco_pensions = . if eco_pensions > 7

rename p2909 eco_environment
replace eco_environment = . if eco_environment > 4

rename p2910 eco_necessity
replace eco_necessity = . if eco_necessity > 4

gen wave = 1

save "Panama_wave1.dta", replace


**************************************

* WAVE 2 Panama
* 2004-2009

**************************************

clear all
import spss using "BASEDATOS_PANAMA_53.sav"
keep nestu pais partido legis p58 p59 p60 p60a p62 p63 p67 p65 p28 p64 p34 p2901 p2902 p2903 p2904 p2905 p2906 p2907 p2908 p2909 p2910

elabel variable (*) ("")

label variable nestu "Study number"
label variable pais "Country name"
label variable partido "Party name"
label variable legis "Legislature"

label drop labels0
label values pais
tostring pais, replace
replace pais="Panama"

replace legis=2009

gen party = "."
replace party = "PN_PRD" if partido ==1 
replace party = "PN_PA" if partido ==2
replace party = "PN_PS" if partido ==3
replace party = "PN_Other" if partido ==4

drop partido 

rename party partido


label variable p58 "Ideology"
rename p58 ID1
replace ID1 = . if ID1 >10
label variable p59 "Ideology of your party"
rename p59 ID2
replace ID2 = . if ID2 > 10

label variable p62 "Female (dummy 1 if woman)"
rename p62 female
recode female (1=0)(2=1)
label define female_label 1 "female" 0 "male"
label value female female_label

rename p60 religious
recode religious (2=0)
replace religious=. if religious > 2
label variable religious "Religious (dummy 1 if believer)"
label define religious_label 1 "believer" 0 "nonbeliever"
label value religious religious_label

gen rel_evangelical=1 if p60a==3
replace rel_evangelical=0 if rel_evangelical!=1
label variable rel_evangelical "Evangelical (dummy 1 if evangelical)"
label define rel_evan_lab 1 "Evangelical" 0 "Other"
label value rel_evangelical rel_evan_lab

gen rel_catholic=1 if p60a==1
replace rel_catholic=0 if rel_catholic!=1
label variable rel_catholic "Catholic (dummy 1 if catholic)"
label define rel_cath_lab 1 "Catholic" 0 "Other"
label value rel_catholic rel_cath_lab

gen rel_other=1 if p60a!=1 & p60a!=3
replace rel_other=0 if rel_other!=1
label variable rel_other "Religious Other (dummy 1 if other)"
label define rel_other_lab 1 "Other" 0 "Evangelical/Catholic"
label value rel_other rel_other_lab

drop p60a

* no question about church attendance
* rename p60b religiosity

label variable p63 "Age (years)"
rename p63 age
replace age = . if age == 99

rename p67 education
label variable education "Highest level of education"
replace education = . if education > 6

rename p65 val_abortion
label variable val_abortion "Opinion about abortion"
replace val_abortion = . if val_abortion > 10

rename p64 val_divorce
label variable val_divorce "Opinion about divorce"
replace val_divorce = . if val_divorce > 10

rename p28 eco_regulated
label variable eco_regulated "Regulate Economy"
label drop labels32
replace eco_regulated = . if eco_regulated > 5
replace eco_regulated = eco_regulated*2
* attention to the rescaling

rename p34 eco_tax
label variable eco_tax "Whether legislator prefers direct or indirect taxes"
label define tax_lab 1 "direct" 2 "indirect"
label value eco_tax tax_lab
replace eco_tax = . if eco_tax !=1 & eco_tax !=2                                                             
rename p2901 eco_prices
label variable eco_prices "Price control"

rename p2902 eco_education_p
label variable eco_education_p "Free primary education"
* combines primary and secondary education

rename p2907 eco_education_u 
label variable eco_education_u "Free university education"
replace eco_education_u = . if eco_education_u > 7
recode eco_education_u (1=1) (2=3) (3=5) (4=7)


rename p2903 eco_housing
label variable eco_housing "Subsidize housing"

rename p2904 eco_employment
label variable eco_employment "Employment creation"
replace eco_employment = . if eco_employment > 7
recode eco_employment (1=1) (2=3) (3=5) (4=7)


rename p2906 eco_health
label variable eco_health "Whether the state should provide health services"
replace eco_health = . if eco_health > 7
recode eco_health (1=1) (2=3) (3=5) (4=7)


rename p2908 eco_unemployment
label variable eco_unemployment "Unemployment insurance"
replace eco_unemployment = . if eco_unemployment > 4

rename p2905 eco_pensions
label variable eco_pensions "Pension"
replace eco_pensions = . if eco_pensions > 7
recode eco_pensions (1=1) (2=3) (3=5) (4=7)


rename p2909 eco_environment
replace eco_environment = . if eco_environment > 4

rename p2910 eco_necessity
replace eco_necessity = . if eco_necessity > 4

gen wave = 2

save "Panama_wave2.dta", replace

********************************************

* WAVE 3 Panama
* 2009-2013

clear all
import spss using "BASEDATOS_PANAMA_71.sav"


keep nestu pais legis partido p29 p3001-p3007 p34a p35 p36 p64 p65 p66 p66a p66b p67 p68 p69 p70 p72

elabel variable (*) ("")

label variable nestu "Study number"
label variable pais "Country name"
label variable partido "Party name"
label variable legis "Legislature"

label drop labels0
label values pais
tostring pais, replace
replace pais="Panama"

replace legis=2013

gen party = "."
replace party = "PN_PRD" if partido ==1 
replace party = "PN_CD" if partido ==2
replace party = "PN_PAN" if partido ==3
replace party = "PN_UP" if partido ==4
replace party = "PN_Other" if partido ==5


drop partido 

rename party partido

label variable p64 "Ideology"
rename p64 ID1
replace ID1 = . if ID1 > 10
label variable p65 "Ideology of your party"
rename p65 ID2
replace ID2 = . if ID2 > 10

label variable p67 "Female (dummy 1 if woman)"
rename p67 female
recode female (1=0)(2=1)
label define female_label 1 "female" 0 "male"
label value female female_label

rename p66 religious
recode religious (2=0)
replace religious=. if religious > 2
label variable religious "Religious (dummy 1 if believer)"
label define religious_label 1 "believer" 0 "nonbeliever"
label value religious religious_label

gen rel_evangelical=1 if p66a==4
replace rel_evangelical=0 if rel_evangelical!=1 & rel_evangelical!=.
label variable rel_evangelical "Evangelical (dummy 1 if evangelical)"
label define rel_evan_lab 1 "Evangelical" 0 "Other"
label value rel_evangelical rel_evan_lab

gen rel_catholic=1 if p66a==1
replace rel_catholic=0 if rel_catholic!=1 & rel_catholic!=.
label variable rel_catholic "Catholic (dummy 1 if catholic)"
label define rel_cath_lab 1 "Catholic" 0 "Other"
label value rel_catholic rel_cath_lab

gen rel_other=1 if p66a != 1 & p66a != 4 & p66a < 8
replace rel_other=. if p66a > 7
replace rel_other=0 if rel_other!=1 & rel_other!=.
label variable rel_other "Religious Other (dummy 1 if other)"
label define rel_other_lab 1 "Other" 0 "Evangelical/Catholic"
label value rel_other rel_other_lab

drop p66a

rename p66b church_attend
label variable church_attend "Attendance to the church (5-point scale, 5 is highest)"
replace church_attend = . if church_attend > 5


* no question about church attendance
* rename p6bb religiosity

label variable p68 "Age (years)"
rename p68 age
replace age = . if age == 99

rename p72 education
label variable education "Highest level of education"
replace education = . if education > 6

rename p70 val_abortion
label variable val_abortion "Opinion about abortion"
replace val_abortion = . if val_abortion > 10

* Attention: how to rescale?
rename p69 val_ssm
label variable val_ssm "Opinion about SSM"
replace val_ssm = . if val_ssm > 2

rename p29 eco_regulated
label variable eco_regulated "Regulate Economy"
replace eco_regulated = . if eco_regulated > 10


rename p34a eco_tax
replace eco_tax=. if eco_tax > 2
label variable eco_tax "Whether legislator prefers direct or indirect taxes"
label define tax_lab 1 "direct" 2 "indirect"
label value eco_tax tax_lab                                                             

rename p3001 eco_prices
label variable eco_prices "Price control"
replace eco_prices = . if eco_prices > 4

rename p3002 eco_education_p
label variable eco_education_p "Free primary education"
replace eco_education_p = . if eco_education_p > 4
* combines primary and secondary education

rename p3007 eco_education_u 
label variable eco_education_u "Free university education"
replace eco_education_u = . if eco_education_u > 7
recode eco_education_u (1=1) (2=3) (3=5) (4=7)


rename p3003 eco_housing
label variable eco_housing "Subsidize housing"
replace eco_housing = . if eco_housing > 4

rename p3004 eco_employment
label variable eco_employment "Employment creation"
replace eco_employment = . if eco_employment > 7
recode eco_employment (1=1) (2=3) (3=5) (4=7)


rename p3006 eco_health
label variable eco_health "Whether the state should provide health services"
replace eco_health = . if eco_health > 7
recode eco_health (1=1) (2=3) (3=5) (4=7)


rename p3005 eco_pensions
label variable eco_pensions "Pension"
replace eco_pensions = . if eco_pensions > 7
recode eco_pensions (1=1) (2=3) (3=5) (4=7)

* Attention to scale
rename p36 eco_natural
label variable eco_natural "Natural resources"
replace eco_natural = . if eco_natural > 10
replace eco_natural = round(eco_natural/2)

gen wave = 3

save "Panama_wave3.dta", replace


********************************************

* WAVE 4 Panama
* 
* Adapted from Valentina

*****************************************************************
{
clear
*Calling the data 2014-2019
import spss using "BASEDATOS_PANAMA_94.sav", clear 
 
*****************************************************************
* Recoding Variables
***************************************************************** 
{
 /*
 

 */

*drop all labels
elabel variable (*) ("")

drop legis 
gen legis = 2019

label drop labels0
label values pais
tostring pais, replace
replace pais="Panama"

* Ideology 
replace ID1=. if  ID1>11
* Ideology Party
replace ID2=. if  ID2>11


// Attendance to church
gen church_attend=RE1b
replace church_attend=. if church_attend>5

gen regular=1 if RE1b==3 | RE1b==4 | RE1b==5
replace regular=0 if RE1b==1 | RE1b==2 

*--------------------------
* Party
*--------------------------
* Only the following parties occur in the data
gen party = "."
replace party = "PN_PRD" if partido ==310
replace party = "PN_MOLIRENA" if partido ==315
replace party = "PN_PP" if partido ==316
replace party = "PN_PANAMEÑISTA" if partido ==319
replace party = "PN_CAMBIO DEMOCRÁTICO" if partido ==320
replace party = "PN_Other" if partido > 339

drop partido 

gen partido=party

*--------------------------
* Abortion
*--------------------------
replace VAL2 =. if VAL2>10
gen VAL2_2=VAL2
gen val_abortion=round(VAL2_2)
lab define VAL221 ///
	1 "Totally against" ///
	10 "Totally agree"
lab val val_abortion VAL221

*--------------------------
* Immigration (not in PANAMA)
*--------------------------
*replace INM101 =. if INM101>6
*gen val_imm=INM101
*lab define VAL21 ///
*	1 "Totally against" ///
*	5 "Totally agree"
*lab val val_imm VAL21

*--------------------------
* SSM
*--------------------------
replace VAL1 =. if VAL1>10
gen val_ssm=VAL1
lab define VAL1 ///
	1 "Totally against" ///
	10 "Totally agree"
lab val val_ssm VAL1

*--------------------------
* Drugs
*--------------------------
replace VAL3 =. if VAL3>10
gen val_drugs=VAL3

lab val val_drugs VAL1

*--------------------------
* Religious
*--------------------------
* Religious label  Are you a  believer?
* Religious
gen religious=1 if RE1==1
replace religious=0 if RE1==2
replace religious=. if RE1>3

lab define bel ///
	0 "Nonbeliever" ///
	1 "Believer"
lab val religious bel

*--------------------------
* Religious Evangelical
*--------------------------
gen rel_evangelical=1 if RE1a==4
replace rel_evangelical=0 if religious==1 & rel_evangelical==1

*--------------------------
* Religious Catholic
*--------------------------
gen rel_catholic=1 if RE1a==1
replace rel_catholic=0 if religious==1 & rel_catholic==1

*--------------------------
* Religious Other
*--------------------------
gen rel_other=.
replace rel_other=1 if religious==1 & rel_evangelical!=1 & rel_catholic!=1


*--------------------------
* Education
*--------------------------
gen education= SOCD7 if  SOCD7<9


lab define uni312 ///
	1 "No education" ///
	6 "Graduate studies"
lab val education uni312

*--------------------------
* Gender
*--------------------------
gen female=1 if SOCD4==2
replace female=0 if SOCD4==1

* Sex label
lab define UNION ///
	0 "Men" ///
	1 "Women"
lab val female UNION

*--------------------------
* Age
*--------------------------
gen age=SOCD5 if SOCD5<100

*-----------------------
* Economy regulate
*-----------------------
gen eco_regulated=EM1
lab define EMI_2 ///
	1 "State" ///
	10 "Market"
lab val eco_regulated EMI_2

*-----------------------
* Free education university
*-----------------------
gen eco_education_u=ROES107 if  ROES107<8
lab define eco_education_u ///
	1 "Against" ///
	7 "In Favor"
lab val eco_education_u eco_education_u

*-----------------------
* State should reduce inequality
*-----------------------
gen eco_inequaltiy=ROES104  if  ROES104<8

lab val eco_inequaltiy eco_education_u

*-----------------------
* State should create employment
*-----------------------
gen eco_employment=ROES103   if  ROES103 <8

lab val eco_employment eco_education_u

*-----------------------
* State companies
*-----------------------
gen eco_companies=ROES101    if  ROES101  <8

lab val eco_companies eco_education_u

*-----------------------
* State wellbeing
*-----------------------
gen eco_wellbeing=ROES102     if  ROES102   <8

lab val eco_wellbeing eco_education_u

*-----------------------
* State health
*-----------------------
gen eco_health=ROES106      if  ROES106    <8

lab val eco_health eco_education_u

*-----------------------
* Pensions
*-----------------------
gen eco_pensions=ROES105       if  ROES105     <8

lab val eco_pensions eco_education_u

*---------------------
* Tax
*---------------------
gen eco_tax = IMP1a if IMP1a < 4
replace eco_tax = 3 if eco_tax > 2

}

// Labeling variables
{
lab var nestu "Study number"
lab var pais "Country name"
lab var partido "Party name (alphanumeric)"
lab var legis "Legislature"
lab var ID1 "Ideology"
lab var ID2 "Ideology of your party"
lab var church_attend "Attendance to the church (5-point scale, 5 is highest)"
lab var religious "Religious (dummy 1 if believer)"
lab var rel_evangelical "Evangelical (dummy 1 if evangelical, 0 if religious but not evangelical)"
lab var rel_catholic "Catholics (dummy 1 if catholic, 0 if religious but not catholic)"
lab var rel_other "Religious Other (dummy 1 if other, 0 if religious but not other)"
lab var val_abortion "Opinion about abortion (1-10 scale, 10 is most in favor)"
lab var val_ssm "Opinion about SSM (1-10 scale, 10 most in favor)"
lab var val_drugs "Opinion about drug legalization (1-10 scale, 10 most in favor)"
*lab var val_imm "Immigrants compete for natives' jobs (1 disagree, 5 agree)"
lab var eco_regulated "Regulate Economy (1-10 scale, 10 agree)"
lab var education "Highest level of education (6-point scale, 1 no education, 6 graduate studies)"
lab var female "Female (dummy 1 if woman)"
lab var age "Age (years)"
lab var eco_inequaltiy "Regulate Inequality between rich and poor (1-7, 7 in favor)"
lab var eco_education_u "Free university education (1-7, 7 in favor)"
lab var eco_employment "Employment creation (1-7, 7 in favor)"
lab var eco_pensions "Pension"
lab var eco_health "Whether the state should own provide health services  (1, 7 in favor)"
lab var eco_wellbeing "Whether the state should guarantee basic wellbeign  (1, 7 in favor)"
lab var eco_companies "Whether the state should own companies  (1, 5)"
lab var eco_tax "Whether legislator prefers direct or indirect taxes (1=direct, 2=indirect, 3=no more tax)"

}
*****************************************************************
* Keep Variables
***************************************************************** 
{ 
	keep nestu pais partido legis ID1 ID2 church_attend religious rel_evangelical rel_catholic rel_other val_abortion val_ssm val_drugs eco_regulated education female age eco_inequaltiy eco_education_u eco_employment eco_pensions eco_health eco_wellbeing eco_companies eco_tax

}

gen wave = 4

*****************************************************************
* Save
***************************************************************** 
save "Panama_wave4.dta", replace
}
*****************************************************************
* Append
***************************************************************** 
{
clear all
append using "Panama_wave1.dta"
append using "Panama_wave2.dta"
append using "Panama_wave3.dta"
append using "Panama_wave4.dta"

save "Panama.dta", replace
*****************************************************************
}


